CALT-68-2812 



Refined matrix models from BPS counting 







Piotr Sulkowsk 



California Institute of Technology, Pasadena, CA 91125, USA 



Abstract 

We construct a free fermion and matrix model representation of refined BPS generating 
functions of D2 and DO branes bound to a single D6 brane, in a class of toric manifolds 
without compact four-cycles. In appropriate limit we obtain a matrix model representation 
of refined topological string amplitudes. We consider a few explicit examples which include 
a matrix model for the refined resolved conifold, or equivalently five- dimensional U (1) gauge 
theory, as well as a matrix representation of the refined MacMahon function. Matrix models 
which we construct have ordinary unitary measure, while their potentials are modified to 
incorporate the effect of the refinement. 
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1 Introduction 



The purpose of this paper is to provide a free fermion and matrix model represen- 
tation of refined topological string amplitudes, and more generally refined BPS counting 
functions in a system of D2 and DO branes bound to a single D6-brane, in toric Calabi- 
Yau manifolds without compact four-cycles. Such a putative free fermion representation 
is interesting, as it would extend earlier results on wall-crossing in D6-D2-D0 system to 
the refined case. The motivation for finding a matrix model representation is as follows. 
In the non-refined case connections between such systems and matrix models are known 
from several perspectives. General relations between topological strings, gauge theories 
and matrix models were postulated by Dijkgraaf and Vafa [1], and related to M = 2 theo- 
ries in [2]. Chern-Simons matrix model for the conifold and generalizations to lens spaces 
were considered in [31 S]. Explicit representation of partition functions of gauge theories 
and topological string theories on corresponding Calabi-Yau manifolds have been found 
in O El [7]. Matrix model representation of partition functions on general toric manifolds 
has been found in [HI |9]. Matrix models encoding wall-crossing phenomena for a class of 
toric manifolds without compact four-cycles have been constructed in [TUl [TT| [T^ . While 
partition functions of four-dimensional gauge theories can be encoded in hermitian ma- 
trix models, a generalization to five-dimensional theories, and more generally topological 
strings on toric manifolds, amounts to considering unitary matrix models [H O El El ITO] . 
All these relations gained new interest with the formulation of the general matrix model 
solution in terms of the topological recursion [13j, and the related remodeling conjecture 
postulated in the context of topological string theory [H]. One might therefore wonder 
if the relation between matrix models and topological strings, and more generally BPS 
counting, extends to the refined well. We also stress that worldsheet definition of 

the refined topological string theory is still not well understood, and the hope that matrix 
model reformulation might give some hint in this context is also an important motivation 
for this work. 

Yet another motivation to study refinement from matrix model perspective arises from 
the AGT conjecture [12]: as proposed in [TB], partition functions of four- dimensional, 
M = 2 theories can be encoded in so-called beta-deformed, hermitian matrix models. Cer- 
tain aspects of this statement were tested in[I7l[l8llia|2ni[2ll[22l[23l[2l[25l|26l[2^ In 
particular the appearance of the beta-deformed measure from the Nekarsov partition func- 
tions has been demonstrated also for both four- and five-dimensional gauge theories and 
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certain topological string theories in |28] . however only to the leading order. On the other 
hand, the formalism of the topological recursion for hermitian models has been extended to 
the beta-deformed case [29j. Therefore one might hope that the refined topological string 
theories could be encoded in unitary, beta-deformed matrix models. However, as explained 
and demonstrated explicitly in ^SU\ \3T\ . this turns out not to be true even in the simple 
example of the resolved conifold. Nonetheless, due to deep consequences of the topologi- 
cal recursion [13], finding some matrix model representation of refined partition functions 
would be quite desirable; such matrix models would presumably arise as some deformation 
of a certain class of already known unitary matrix models. This is the task we cope with 
in this paper, not only from the viewpoint of topological string amplitudes, but also more 
generally in the context of BPS counting and wall-crossing phenomena. The refined matrix 
models which we find involve matrices of infinite size and have ordinary, unitary measure, 
while their potentials are modified in a way which encodes the refinement. We stress this 
is opposite to the beta-deformed models, whose measure is modified, however potentials 
are the same both in refined and non-refined cases. One immediate advantage of our result 
is the fact, that the topological recursion for models with undeformed measure [131 [H] 
is much simpler and tractable than in the beta-deformed case [2^], and could be readily 
applied to gain more insight into properties of refined amplitudes. 

We recall that there are various definitions of refinement whose physical equivalence is 
not quite clear, however the agreement of the resulting exact solutions is a strong argument 
for an underlying common, general structure. In all these so-called refined theories a 
dependence on a single parameter, such as string coupling Qs or the background h in gauge 
theories, is replaced by a dependence on two parameters, customarily denoted ei and €2- 
In the context of gauge theory refined amplitudes arose from their formulation in the Q- 
background [32] • In case of topological strings on non-compact, toric manifolds, refinement 
was introduced in terms of refined BPS counting, reformulated combinatorially in terms 
of the refined topological vertex |33[ and shown to agree with gauge theory results 
in Q background in [35| [36] . From the viewpoint of AGT conjecture refined amplitudes 
are encoded in relevant conformal blocks of two-dimensional CFT, and the corresponding 
beta-deformed matrix models are characterized by the Vandermonde determinant raised to 
the power /3 = — ei/e2. In the context of wall-crossing and BPS counting in a system of D6- 
D2-D0 branes on toric manifolds, following and in parallel with non-refined developments 
in [371 l38| [39| H0| WT\ B2| B3| HI], refined amplitudes were considered from physical and 
mathematical perspectives in [15| Among multitude chambers in which (refined) 
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generating functions of D6-D2-D0 bound states are known, there is a special chamber in 
which they agree with topological string amplitudes on the same Calabi-Yau manifold, and 
in particular the agreement with the refined topological vertex calculation was shown in 
|46| . This is this last formulation of the refinement on which our derivation is based. 

To find refined matrix models we follow a strategy which extends a non-refined pre- 
sentation of [lOjfl Firstly, generalizing the results of |43j, we construct free fermion repre- 
sentation of crystals representing refined BPS states in question. This allows to write the 
refined BPS generating functions Z^^^ in a chamber specified by n as 

z::^ = {^':^\w:^\^-f), (i) 

where \^^±^) are states representing a manifold in question, and W'^^^ are wall-crossing 
operators which determine a chamber of interest. Then we turn these fermionic correlators 
into a unitary matrix model form. The refined character of fermionic correlators results 
in a modified form of matrix model potentials. Similarly as in O El [10], our potentials 
have nontrivial string coupling dependence to all orders. While our results are valid in all 
chambers, in the so called commutative chamber we obtain matrix model representation 
of refined topological string amplitudes. 

To briefly exemplify our results, we recall first that the refined topological string ampli- 
tude for the resolved conifold with Kahler parameter Q (or equivalently five-dimensional, 
U{1) gauge theory) is given by 

oo 

zr4 = M{t,M)X{{i-Qt't'%), (2) 

fc,Z=0 

where ti = e-^\t2 = e^^, and M{ti,t2) = llZ=o^^ " ^i^^4)"^ is the refined MacMahon 
function. To find a matrix model representation of Z^^p > ^^st construct a general refined 
BPS generating function in the form ([1]), where in the case of the conifold n is a single 
integer. We then translate such a fermionic correlator into a matrix model form, and in 
— )■ oo limit, which corresponds to the so-called commutative chamber, we find the matrix 
model representation (written in terms of eigenvalues = 6*"*=) of the refined topological 
string amplitude 

-^ref _ f ^ttTJ TT + ^fc^l'^ ) + 4/ ^fe) 



^Our results were obtained independently and before an overlapping work |47l I48| appeared. 



4 



where VU is the ordinary unitary measure (see (1181) ). To the leading order the above 
integrand gives rise to the following potential 

Viu; /3) = ^u' - (1 - /3-i)Li2(-e^") - Li2(-ge-^") + Oig„ (3) 

In what follows we also present matrix models associated to other chambers of the Kahler 
moduli space. We can also immediately note that in the limit Q — )■ 0, the above result 
reduces to a matrix model representation of the refined MacMahon function, with the 
exact integrand given by a deformed theta-function, which in the genus expansion gives 
a /3-deformation of the gaussian potential of the Chern-Simons matrix model (such that 
both the dilogarithm term, as well as 0{gs, (3) corrections, vanish for /3 = 1). In the main 
text we discuss in more detail other explicit results for C'^, conifold, or resolution of C'^/Z2 
singularity. Similarly as in [12j we postulate a relation of those refined matrix integrands 
to open BPS amplitudes. 

The paper is organized as follows. In section [2] we recall definitions and basic properties 
of refined BPS invariants and introduce relevant notation. In section [3] we extend formal- 
ism of P3] to the refined case and present fermionic representation of refined generating 
functions. In section H] we turn these refined fermionic results into matrix models and 
describe their properties. Section [5] contains a discussion. 

2 Refined wall-crossing in D6-D2-D0 system 

Refined degeneracies of D2 and DO-branes bound to a D6 brane on a Calabi-Yau man- 
ifold X can be encoded in a generating function 

a, 7 

with DO-brane charge represented by a G Z, D2-brane charge represented by 7 G H2{X, Z), 
and a chamber in the Kahler moduli space specified by (possibly a set of parameters) n. 
Let T-La^-yin) denote a space of BPS states with given charges 0,7 and asymptotic values 
of moduli corresponding to a chamber n, and J3 denote a generator of the spatial rotation 
group. For fixed charges a, 7 and a choice of chamber n, refined degeneracies 

f^S(n;2/) = TrK„,.(n)(-l/)''^ (4) 

are interesting invariants if X does not posses complex structure deformations, which is the 
case for non-compact, toric manifolds which we consider in this paper. These invariants 
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were argued in P3] to agree with motivic Donaldson-Thomas invariants of |49] . and in 
the case of the resolved conifold the corresponding BPS generating functions were derived 
using the refined wall-crossing formula, and encoded in a refined crystal model. From 
mathematical viewpoint, and in terms of dimer models, such analysis was extended to quite 
a general class of toric manifolds without compact four-cycles in |3B], and shown therein 
to agree, in the commutative chamber, with refined topological vertex computations. For 
?/ = 1 all these invariants reduce to ordinary non-refined invariants, whose generating 
functions were encoded in dimer or crystal models in [39| H0| HI], and represented in 
the free fermion formalism in HI]. In the next section, based on definitions of BPS 
generating functions in terms of dimers or crystals constructed in I46 | . we will extend 
such free fermion formalism to the refined models. 

Before proceeding we present in more detail a class of manifolds we are interested 
in. Similarly as in P5| [TU], we consider toric, non-compact Calabi-Yau manifolds without 
compact four-cycles, whose toric diagrams arise from a triangulation of a strip. Such a 
diagram consists of + 1 vertices, and there are P^'s in the geometry with Kahler 
parameters denoted Qp = e~"^p, p = 1, . . . , N. To each vertex in the diagram we associate 
a type Tj = ±1. If the local neighborhood of P^, represented by an interval between vertices 
i and i + 1, is C(-2) © O, then r^+i = r^; if this neighborhood is of C(-l) © type, 
then Tj+i = — Tj. We choose a type of the first vertex as ri = +1. 

We also need to introduce relevant notation for refined quantities. In the non-refined 
case the string coupling is related to the DO-brane charge as g = e"^^ The refinement 
is encoded in an additional parameter /3. Instead Qs and /3 it is more convenient to use a 
pair of parameters 



The variable ?/ in (jl]) can be expressed as y = ti/q = q/t2, so that y"^ = ti/t2 = ■ In this 
notation the non-refined limit y = 1 corresponds to /3 = 1, for which ei = —€2 = Qs and 
ti = t2 = q and B = 0. 



so that /? = ——, 6162 = —g[ 



g. We often use the exponentiated counterparts 



ti = e 



and also introduce 




6 



With the above notation we can present some exphcit BPS generating functions whose 
matrix model representation we are going to find. The simplest manifold one can consider 
is C^, for which one gets the refined MacMahon function [33j . see fig. [H 

oo ^ 

Z'^' = M{t„t,)=l[--^. (5) 

k,l=0 ^ '^1 

In this case there is no Kahler parameter, and therefore there are no interesting wall- 
crossing phenomena. 

We note that one could consider more general family of refinements parametrized by 
S, such that Ms{ti,t2) = nfc°/=o ~ ^i'''^^ ^4 ^ ) ^- For simplicity, in what follows we 
choose the value 5 = 1 (note that in 05] another choice 6 = was made). 

The resolved conifold provides a basic non-trivial example of wall-crossing, with a 
set of chambers parametrized by an integer n (in the refined case one might also consider 
additional invisible walls, which we do not discuss here). Corresponding refined generating 
functions were computed in [l5] using a refined wall-crossing formula and in the chamber 
labeled by n — 1 they read 

oo 

zzi^"'' = M{t,M)\ n (1 - ^^^'4)) ( n (1 - (6) 

k,l=0 k>l, l>0, k+l>n 

In the commutative chamber n — oo the terms in the last bracket do not contribute 
anymore and the BPS generating function is simply related to the refined topological 
string amplitude given in ([2]) 

On the other hand, in the non- commutative chamber n = 0, the refined generating function 
is given by the modulus square of the refined topological string amplitude. 

For a resolution of C^/1j2 singularity there is also a discrete set of chambers 
parametrized by an integer n and the corresponding BPS generating functions read 

oo 

zff^=M(t„t,)^(n (i-Qtr4)")( n {i-Q^'t%)-'). (7) 

k,l=0 k>l, l>0, k+l>n 

It is harder to write down generating functions for arbitrary chamber of an arbitrary 
geometry of our interest. However this can be done in for the non-commutative chamber of 
arbitrary geometry, where - similarly as in the non-refined case - BPS generating function 
is given by the modulus square of the refined topological string amplitude 

Zq'^^ = l^top^l^ = 2^top {Qi)2t^ {Qi ^). (8) 
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The (instanton part of the) refined topological string amplitude is given by [551 [55] 

oo 

z::f{Q,)=M{h,t,)^ n n {^-iQ^Q^+l■■■Q^--l)tl^%y"^\ (9) 

k,l=0 l<i<j<N+l 

with the notation introduced above. 

3 Refined wall-crossing and free fermions 

The problem of counting bound states of D6-D2-D0 branes for local toric Calabi-Yau 
manifolds without compact four-cycles has been formulated in the free fermion formalism 
in 1^ IS]. Among many advantages of such a representation is its immediate relation to 
melting crystals, as well as to matrix models, which was exploited in [TOl [12]. Here we 
wish to extend such free fermion formalism to capture refined BPS invariants, as defined 
in [13 [ig. 

We consider first statistical models of colored pyramids. In the non-refined case [13], to 
a geometry consisting of P^'s one associates a crystal which is sliced into layers in A^-|- 1 
colors, denoted qo, qi, q2, ■ ■ ■ , qN- In the non-refined case, parameters qi, . . . ,qN encode 
Kahler parameters of the geometry Qi, . . . ,Qn, while the product j^^^o'?* mapped to 
(possibly inverse of) q = e~^\ In the refined case the assignment of colors is more subtle, 
as it must take into account a refinement of a single parameter q into ti and t2 introduced 
above. In particular, in the non-commutative chamber gj^o are mapped (up to a sign, as in 
the non-refined case) to Qi, however we will have to replace go by two refined colors ^q^^ or 
ql^\ so that tj = qQ^i ■ ■ ■ qn^ for ^ = 1? 2. The simplest case of refined plane partitions, 
discussed also in [33], is shown in fig. [1] For other manifolds, in other chambers we will 
find more complicated assignment of colors. 

In [13] the structure and coloring of a given crystal, corresponding to a particular toric 
geometry, was encoded in fermionic states so that the generating function of BPS 

invariants could be written as a superposition of two such states (with additional insertion 
of wall-crossing operators in non-trivial chambers). In this section we construct a refined 
states with similar properties. In the non-commutative chamber the states which 

we construct are such that 

zi^'f = (10) 

re/ 

We also construct refined version of wall-crossing operators , such that the BPS 
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Figure 1: Refined plane partitions which count D6-D0 hound states in C^, as 
seen from the bottom (i.e. a negative direction of z- axis). Stones in each layer 
which intersects a dashed or solid line have weight ti or t2 respectively. The 
resulting generating function is the refined MacMahon function M{ti,t2). 



generating function in n'th chamber can be written as 

^ref ^ ^^f\wl"^\n^[^f). (11) 

In section [XT] below we construct states \^±'^) for arbitrary manifold in a class of our 
interest. In section [3l2] we construct states \^±'^) and wall-crossing operators W^^^ for 
all chambers of the resolved conifold and a resolution of C^/Z2 singularity. We follow 
conventions used in [121 [IDl [12] , which are summarized for convenience in appendix [XI 

3.1 Arbitrary geometry — non-commutative chamber 

In this section we construct fermionic states \^±'^), which allow to write the BPS 
generating functions in the non-commutative chamber as claimed in ([TOj) . Similarly as 
in the non- refined case, the states are constructed from an interlacing pattern of 

vertex and weight operators. As the refinement does not modify the three-dimensional 
shape of the corresponding crystal, the assignment of vertex operators is the same as in the 
non-refined case jl3] and can be similarly read off from the toric diagram. In particular, to 
the i'th vertex in the toric diagram (of type Tj given above) we associate a vertex operator 
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r^(x), such that 

Examples of this assignment for C'^, conifold, and a resolution ofC^/Z2 singularity are 
shown in fig. [2J 




Figure 2: Toric diagrams and assignment of vertex operators in case of 
(left), conifold (middle), and a resolution o/C^/Z2 singularity (right). A sign 
(B or Q in each vertex denotes a corresponding type Xj = ±1. 



The structure which is modified in the refined case is the assignment of colors, which are 
encoded in the weight operators. A product of such operators T'^{x) is interlaced with 
weight operators in the following way. We introduce operators Qi representing colors 
Qi, for i = 1, . . . ,N, and in addition two other colors q^^^ and which are eigenvalues 
of Qo^^ and Q'^^ ■ Operators Qi, ... ,(5 at are associated to in the toric diagram, and we 
define 

Q^'^ = &o^Qi---Qn, U = qfqi---qN, for 2 = 1,2. (12) 
Now we introduce 

= T^^{x)Q,V^_^{x)Q2---T^^{x)QnT^^^\x)Q^^\ 
A.{x) = n{x)QiV^^{x)Q2---r^{x)QMT^_!'+\x)Qf. 

Commuting all Qj's to the left or right we also introduce 

A,{x) = &'TA,{x) = r^[xt,)T^AT^^{^^ 

A^x) = A4x) (g('))"' = ^L^a:)^L2(a;gl)^!:'(xglg2)■■■^?^+'Nlg2•••gJv)• 
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When the argument of any of these operators is x = 1 we often use a simphfied notation 

Finally we can associate to a given toric manifold two states 

{ni^^l = {0\...A4l)A4l)A4l) = {0\...A4tl)A4tMM, 

= A4i)A_(i)A4i)...|o) = A_(i)A_(t2M_(t^)...|o), 

where |0) is the fermionic Fock vacuum. 

Our first claim is that the refined BPS generating function can be written as 

z;^^ = {ni'f\n^jf) ^ z,,,m2top{Q-'), (is) 

with ZtopiQi) given in and under the following identification between g-j parameters 
which enter a definition of \^±'^) and string parameters Qi = e~^' (for 2 = 1,..., A^) 

Qi = {TiTi+i)Qi, 

and with refined parameters ti^2 identified as in (fT2|) . This result, in the special case of C^, 
conifold and C'^/Z2 geometries, reproduces formulas ([5]), and ([7]). 

To prove (IT^ for general geometry, we first note that commuting operators A^{x) with 

My) 

A+{x)A_{y) = A_{y)A+{x) C{x,y), 

gives rise to a factor 

C{x,y)= . .^+1 n ((l-(r.r,)a;yti(gigi+i...g,_i))(l--|^^^^^^)) 

Now we write the states in terms of A± operators, and commute r± within each 

pair of and A_ separately 

oo oo oo 

Zo'' = = (oi(n^+(ti))(n^-(^2))io) = n c^^^- 

i=0 j=0 i,j=0 

This last product reproduces modulus square of the refined topological string partition 
function in (fT3|) and therefore proves the claim (fTOj) . 
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3.2 Conifold and C^/Z2 all chambers 

Fermionic representation can also be extended to non-trivial chambers. Even though 
this can be done for general geometry without compact four-cycles, for simplicity we restrict 
our considerations to the case of a conifold and a resolution of C'^/Z2 singularity, which 
involve just one Kahler parameter Qi = Q. In those cases, in a chamber labeled by — 1 
we find the following representation of BPS generating function 

zZ{ = {^7'\w':!,\nL^'), (14) 

ref 

where W^_i represents appropriate wall-crossing operator. In these both cases the toric 
diagram has two vertices, the first one of type ri = 1 and the second one denoted now 
T = r2, and r = =1=1 respectively for the conifold and C'^/Z2. A crystal associated to 
the expression fll4p has n stones in the top row and can be sliced into interlacing single- 
colored layers. The assignment of colors is analogous as in the pyramid model discussed in 
[15| HH] (however our convention is slightly different, and corresponds to integer and non- 
symmetric, rather than half-integer and symmetric powers of ti 2 in [45j). The pyramid 
crystal for the conifold is shown in fig. [31 The coloring and weights for C^/Z2 are the same 
as for the conifold, even though the plane-partition shape of C^/Z2 crystal is different 
(though very analogous) than pyramid-like conifold crystal, see fig. |H Using the notation 
introduced above, the assignment of colors is determined as follows. 
All stones on one side of the crystal are encoded in the state 

= (o| . . . (r+(i)gir;(i)go) (r+(i)gir;(i)go). 

The Kahler parameter g, as well as the parameter ti, are determined as 
qi = ^gtj"", qo = T-^, so that go^i = ^i- 

Then the extended crystal, which has n — 1 additional stones in the top row, is con- 
structed by an insertion of the operator 

= (r_(i)gir;(i)goP)(r-(i)Qi?r;(i)gof^)--- 
... (r4i)gig^^r;(i)gog(^^)- 

This operator consists oi n — 1 terms of the form ^r_(l)gig*^r^(l)goq'^(*+^)^ j for i = 
0, . . . ,n — 2, where in each consecutive dark or light slice of stones we insert one additional 
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Figure 3: Refined pyramid crystal for the conifold, in the chamber corresponding 
to n stones in the top row. Along each slice (as indicated by broken or solid 
lines) all stones have the same color, assigned as follows. On the left side 
(along broken lines), each light (yellow) and dark (red) slice has color denoted 
go and qi respectively. Moving to the right, in the intermediate region (along 
solid lines), a color of each new light or dark slice is modified by respectively 
q^^ factor (with respect to the previous light or dark slice). On the right side 
(again along broken lines), each light or dark slice has again the same color, 
respectively qoQ^^^ or qiq^^^^^\ The assignment of colors in the intermediate 
region (along solid lines) interpolates between constant assignments on the left 
and right side of the pyramid. 



operator q^^, which changes the weight of each stone in this shce by q"^^ = (^1/^2)^^ (with 
respect to the previous shce of the same hght or dark color). 

Finally, all stones on the right side of the crystal have again the same light or dark 
color, and the corresponding state reads 

IQ'-y) = (r_(i)Qig(^^r:(i)goF^) (r_(i)gig(^^r:(i)Qof^) • • • |o). 

Therefore the varying weights in the middle range (along solid lines in fig. [3] and H]) 
interpolate between fixed weights of light and dark stones on two external sides of a crystal. 
We can now commute away all weight operators in the above expressions, using relations 
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Figure 4: Refined pyramid crystal for the resolution of C^/Z2 singularity, in 
the chamber corresponding to n stones in the top row, as seen from the bottom 
(i.e. a negative direction of z- axis). Even though the three-dimensional shape 
of the crystal is different than in the conifold case, the assignment of colors is 
the same, see fig. 0. 



from appendix \^ This leads to the representation 

oo n— 2 oo 

= (o| ( n r+(it)r;(tt/gi)) ( n r4t2)r;(gr V)) ( n 

k=l i=0 fc=0 

Finally, commuting all vertex operators, we find 

oo oo 

zi'i,=M{t,Mf n n {^-Q^'t%)-\ (15) 

fc=lJ=0 k>l,l>0,k+l>n 

where r = =pl respectively for the conifold and C'^/Z2. This indeed reproduces ([6]) and 
([7]), and agrees (up to half-integer convention of ^1^2) with the results of [45|. 

4 Matrix models from fermions 

Once the generating function of Donaldson-Thomas invariants is written in the 
fermionic formalism, it can be turned into matrix model form upon inserting appropri- 
ately chosen identity operator in the correlator (HM 

= {n':^\iw:!,\Q'y). (16) 
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The identity operator I is represented by the complete set of states |-R)(-R| (representing 
two-dimensional partitions). We can use orthogonality relations of [/(oo) characters XRi 
and representation of these characters in terms of Schur functions xr = sr{z) for z = 
{zi,Z2,Z3, . . .), to write 

R P,R P,R 

= ['^u[l[T'_{z,m)(^{0\l[T'4z^')), (17) 

k k 

where VU denotes the unitary measure written in terms of eigenvalues 

VU = \[duk\[\zu-z,\\ Zk = e'''K (18) 

k k<j 

The identity operator in the above form can be inserted into (ITB]) which results in an 
expression involving only vertex operators r|^^\ Then we can commute vertex operators 
away, again using relations from appendix |3 which leads to a matrix model with the 
unitary measure VU. In the non-commutative chamber all factors arising from commuting 
these T^^^^ operators depend on Zk and contribute just to the matrix model potentials. In 
other chambers additional factors may arise which are independent of z^, and which, in a 
chamber labeled by n, contribute to some overall factor /„. Thus in general we write the 
Donaldson-Thomas generating function as a matrix model in the form 

= fn I VUHe-'^^^^-^K (19) 

k 

and it is convenient to introduce a factor in front of the potential V{z;l3), or work 
with a rescaled coupling Qsl^'^^"^ ■ 

In fll6p the identity operator has been inserted in a specific location. In fact there is 
large freedom of where this insertion should be chosen, which leads to various form of a 
matrix integrand. In [12] it has been shown that those various integrands can be identified 
with open BPS generating functions in various open chambers. We will also comment on 
a possible similar interpretation of refined integrands in what follows. However, let us 
first restrict to the a specific choice (fT6|) and discuss resulting matrix models. We use the 
following notation for a deformation of a theta-function 

oo 

Q{z-MM) = X{{l + zt{^'){l + ti/z) 

j=0 

to express certain integrands of matrix models that we come across. 
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4.1 Arbitrary geometry non-commutative chamber 

As the first explicit example, we find matrix model representation of the refined BPS 
generating function in the non-commutative chamber. We start with the expression (fT6|) 
with no Wl^li insertion, and use the form of \^'±'^) derived in section 13.11 Performing 
the computation described above we get, in the non-commutative chamber for general 
geometry, the following matrix model: 

r ^ 

i.e. we identify e"^^'-^''^-* = Hilo ® ' ' ' '^2)^'^^- The product over / runs 

over all vertices and in this chamber we identify Kahler parameters Qp with weights Qp via 

% = {'TpTp+l)Qp- 

Some special cases of the above result include: 

• for the generating function Z'^'^^ = M{ti,t2) is given by the refined MacMahon 
function and we find that the corresponding potential is a refined theta function 

00 

e-^v-C.;/?) = Y[{1 + zt{'')il + ti/z) = Q{z; t,M) (20) 
i=o 

• for the conifold the non-commutative generating function Zq°^^^°^'^ determined from 
(E]) gives rise to a matrix model with the following potential term 

• for C'^/Z2 the non-commutative generating function Zq determined from ([7]) gives 
rise to a matrix model with the following potential term 

e-^y^--^) = \{{l + zt{^'){l + ti/z){l^zt{^'/Q){l+'^) = 

j=0 ^ 

= Q{z-tiM)Q{z/Q-hM) 

4.2 matrix model 

Let us consider the simplest refined matrix model, corresponding to geometry, with 
the exact potential given in (l20l) . Similarly as in [HI [50], one might expect that its behavior 
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is governed by the leading order term in the potential Using the asymptotics 
this leading behavior reads 

g-^y(^;/3) _ [_i(log^)2-(l-/3-i)Li2(-.)+0(g„/3)] _ ^21) 

The first, quadratic term in the potential is the same as in the non- refined case. The term 
involving Li2(— 2;), as well as all higher order terms 0{gs,P), vanish for /3 = 1. Therefore, 
for /3 = 1, we obtain a Chern-Simons matrix model which indeed is known to give rise to 
MacMahon function in — )■ 00 limit ^ilOj. For general /3, a resolvent u){p) for a unitary 
model with the above potential can be found using the Migdal integral, as discussed in 
detail in [TU| \5T\ . This requires bringing the measure into a hermitian Vandermonde 
form, which introduces additional Tlogz term in the matrix potential, with the 't Hooft 
parameter T = {gs(3~^^'^)N. For the lowest order terms of the potential arising from (12T|) . 
this leads to 

a,V(z, = + (22) 

Assuming a one-cut solution of the matrix model, and in terms of the rescaled coupling 
dsP'^^"^, the Migdal resolvent is then given bj|^ 



1 f dz dzV{z) \/{p-a-){p 



27ri p- z ^[z - a-){z - a+) 
so that the endpoints of the cut a_ and are encircled counter-clockwise by the integration 
contour. Moreover, one has to impose the following consistency condition on the resolvent 

lim u{p) = -, 

p->-oo p 

which imposes certain conditions on the end-points of the cut a±. We find that for the 
potential (122]) these conditions take form 



e^/^ (23) 



, ( , ^ ■ =) = 2e-^/l 24 

V + + ^(a^ + i)a_ / 



A useful result [10^ in such computations is ^ f d^ ^°e{z+c) V(p >>) 



2T J 2Tvi z(p-z) ^(z-a){z-b) 

- 1 log (V( ^^ppV«'pKa-P) Y V(P^a)iP-b) iog( V(^-V^^ )\ This arises from con- 



tour integrals around poles z = and z = as well as along the branch cut of the logarithm (—00, — c) 

which is found using f ^ = (V(.-a)(b-p)-V(.-b)(a-p))- _ 

(3:-p)V \J {p-a){p-V) (p-x)^ {p-a){p-b) 
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For the non-refined case /3 = 1 these equations simphfy and can be exactly solved [TU]. For 
arbitrary /3 the cut end-points found in [TU] get corrections in (1 — 

a± = -1 + 2e-^ ± 2ze-^/Vl - + C(l - 

which leads to a /3-deformed spectral curve. To find these corrections (9(1 — in the 
exact form appears not easy, and it would be interesting to compare the resulting curve with 
the quantum curve of the beta-deformed formalism of j29]. In particular they both give 
rise to the same result in the four- dimensional limit |I31], so understanding a discrepancy 
of the five-dimensional results is an important issue. It would also be interesting to find 
the partition function for the above model with finite 't Hooft coupling T, and verify if it 
is related to the refined conifold topological string amplitude, as is indeed the case in the 
non-refined case. 

As already mentioned before, we can also obtain more general matrix models by insert- 
ing the identity operator in various places in the fermionic representation of BPS function. 
In particular, inserting it at position A; in a string A_ operators in (fT3|) in case, we get 
the following representation 




with the prefactor = )f:2) Il^fc ni^o(-'- ^1^2)- -'■^ non-refined case in [12j . 

the above integrand with identification ti = t2 = q was related to open BPS generating 
function in an open chamber labeled by k. It would be interesting to extend such an 
interpretation to the refined case too. In particular, we note that in the limit k — > 00, 
which should correspond to the ordinary open topological string amplitude, the above 
integrand indeed reduces to one particular form of a refined brane partition function in C'^ 
computed in [33j. 



4.3 Conifold all chambers 



Using the representation (1161) and fermionic results found in section 13. 2[ we find the 
following matrix model for the conifold in the n'th closed BPS chamber (corresponding to 
a pyramid with {n + 1) stones on top), and with identity representation (??) inserted at 
position /c in a string A_ operators in ( IT3|) 
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with relevant prefactor fn,k{<l, Q)- In the special case = we get 

oo oo 
k=l,l=0 k>l,l>0,k+l>n+l 



with the prefactor given by 



k j=0 



/n(g,Q)=(nnT37^)(n n a-n^/Q) 

i=l A:=0 ^ 2 4=1 j=n+l-i 

We determine cut end-points, which would follow from deformed equations analogous to 
(??) and (??), and which would lead to a /3-deformed mirror curve; however this is quite 
involved technically. 

On the other hand, it is instructive to analyze the limit of commutative chamber n ^ oo 
. Firstly, in this limit the prefactor simplifies to f^o = M{ti,t2). We also know that the 
total value of the BPS generating function reduces in this limit to the refined topological 
string amplitude, up to a single factor of M{ti,t2). In consequence in the commutative 
chamber we get a matrix model representation of the refined topological string amplitude 
for the conifold 



k,l=0 



A 111 {l + tiQ/z,) 

Now the leading order potential term is a modification of the potential fl2ip by a Q- 
dependent dilogarithm term 

Viz; /3) = -^(log zf - (1 - r')Li2(-^) - U2i-Q/z) + 0{gs, /?). (26) 

In the limit Q — the above topological string partition function becomes just the refined 
MacMahon function, and the matrix integral consistently reproduces C'^ result fl20l) . 
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4.4 C^/Z2 all chambers 

The results for C'^/Z2 arise similarly as those for the conifold. Using the representation 
(flGjl and fermionic construction from section 1X2) we find 

oo oo 

= M{uMf n i^-QAAr' n {I - Q-'t\t'^-' = 

fc=l,Z=0 fc>l,/>0,fe+/>n 

/oo 

where 

Tl OO Ti OO 

/«(^.«-(nnT4i)(n,n,rrii7i^)- 

1=1 fe=0 1=1 j=?i+l— I i ^/ ^ 

In particular, in the commutative chamber n — )■ oo we get again /oo = M(ti, ^2)- Therefore 
in the commutative chamber we get a matrix model representation of the refined topological 
string amplitude 

00 ^ „ 00 

k=l,l=0 ^ 1 2 J ^ j^Q 

In the limit Q — > we again recover the refined MacMahon function, as well as the 
expected integrand of C'^ matrix model (120]) . In this case it is also straightforward to find 
more general matrix models, analogous to f l25p . which would presumably be related to 
refined open amplitudes. 



5 Discussion 

In this paper we have found a free fermion, as well as a unitary matrix model rep- 
resentation of refined BPS generating functions of DO and D2-branes bound to a single 
D6-brane, and in particular topological string amplitudes, in toric Calabi-Yau manifolds 
without compact four cycles. We mainly considered explicit examples of C'^, conifold, 
and C^/1j2 geometries, as well as an arbitrary geometry in the non- commutative chamber, 
however generalization to other chambers for manifolds in this class is straightforward. A 
general consequence of our results is the fact that refined generating functions, at least for 
the class of manifolds which we considered, have nice properties of ordinary matrix model 
expressions [T3| [Ti] . such as integrability, symplectic invariance of associated free energy 
coefficients Fg, automatic appearance of the whole family of differentials W^, etc. One 
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advantage of our representation is that these properties are much better understood for 
ordinary matrix models, rather than for matrix models for beta-ensembles [.29J, which in 
fact are known not to reproduce the refined topological string amplitudes [30t [3Tj . It is 
also important to understand a difference between these two beta-deformations. As follows 
from the results of [3T], in case of the conifold (or five-dimensional U{1) gauge theory), the 
four-dimensional gauge theory limits of both deformations agree. Understanding the origin 
of a discrepancy in five- dimensional deformation should lead to interesting new insights. 

There are many other questions which require further investigation. Firstly, a nontrivial 
task is to find spectral curves of our models. As we discussed, these would be /3-deformation 
of curves found in a non-refined case in [10]. Having known such curves would allow to 
apply the topological recursion to recover quantities and Fg explicitly from matrix 
model perspective. This appears nontrivial, in particular due to all order Qs corrections to 
our potentials. However these corrections arise from terms involving quantum dilogarithms. 
Potentials which involve quantum dilogarithms were considered also in 13120], where it was 
shown that higher Qs essentially do not modify resulting invariants, and one can effectively 
consider a leading order contribution to the potential, similarly as in (12T]) and ( 126|) in 
our case. It would be interesting to confirm if analogous phenomenon takes place for the 
potentials which we consider. 

Furthermore, it would be interesting to extend our discussion to the open string case, 
on one hand refining the discussion in [12j and providing M-theory derivation of putative 
open BPS generating functions, and on the other relating W-( to brane amplitudes in 
matrix models in the topological string limit. In particular this should provide a deeper 
understanding of nontrivial prefactors in intermediate chambers. 

It would of course be interesting to extend our results to toric manifolds with compact 
four-cycles, in particular those related by geometric engineering to gauge theories. This 
might be possible by considering more involved crystal models, such as those in [52]. 

Among other questions, it is interesting what our matrix models compute for finite size 
of matrices A^. It was shown in [10] that in the non-refined case finite A^ engineers more 
complicated toric manifolds with an additional two-cycle (as is already the case in the 
Chern-Simons matrix models |1], where a finite 't Hooft coupling encodes the size of the 
single of the resolved conifold). In particular it is tempting to speculate whether the 
matrix model (l20l) with finite A^ would also provide the refined conifold topological string 
partition function. 

It would also be interesting to understand the issues of holomorphic anomaly and 
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modularity and make contact with discussions in [5U| IM] , and more generally with 
large literature on refined invariants. 

We hope that continuing this line of research would be a rewarding experience. 
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A Free fermion formalism 



In this appendix we summarize free fermion formalism used in |^ [T0| [T2] . We consider 
the Heisenberg algebra [am-, «-n] = nSm,n and define vertex operators 



which satisfy commutation relations 

I — xy 1 — xy 

r;(x)r_(y) = (1 + xy)T^{y)T'4x), T+ix)T'_iy) = (1 + xy)T'_iy)T4x). 

These operators act on fermionic states corresponding to two-dimensional partitions 

fj,, as 

r-(^)l/i) = J^xl^l-l'^llA), r+(a;)|/.) = 5^xl-H^I|A), (27) 

r'-WI-") = E rv(x)i^) = J2 (28) 

where -< is the interlacing relation. We also consider various weight operators Qg, with 
eigenvalues representing colors and denoted qg, such that 

4|A) = gl'l|A), 
and their commutation relations with vertex operators read 

T+ix)Qg = QgT+{xqg), T'4x)Qg = QgT'^xqg), (29) 

QgT4x) = T_{xqg)Qg, Q,r'_(a;) = T (xg,)Q,. (30) 
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